Ratchet behavior in nonlinear Klein-Gordon systems with point-like inhomogeneities 
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We investigate the ratchet dynamics of nonlinear Klein-Gordon kinks in a periodic, asymmetric 
lattice of point-like inhomogeneities. We explain the underlying rectification mechanism within a 
collective coordinate framework, which shows that such system behaves as a rocking ratchet for point 
particles. Careful attention is given to the kink width dynamics and its role in the transport. We 
also analyze the robustness of our kink rocking ratchet in the presence of noise. We show that the 
noise activates unidirectional motion in a parameter range where such motion is not observed in the 
noiseless case. This is subsequently corroborated by the collective variable theory. An explanation 
for this new phenomenon! is given. 



I. INTRODUCTION 

The study of transport mechanisms at the mesoscale 
level is of great importance nowadays. Specifically, the 
so-called ratchet systems have shown to be proper can- 
didates for explaining unidirectional motion to biological 
systems Q , and have important physical applications for 
nano- and micro-scale technologies H H ||. Many of 
these models have been developed in the simplepicture 
of point-like particles 0, 0, ( see the reviews @ |j| for 
details) . Such scenario has been subsequently generalized 
to spatially extended systems 0>O>0> where much at- 
tention has been paid to situations where the net motion 
arises through time-symmetry breaking [l3L ITlj . This 
kind of ratchet phenomenon has been recently observed 
in long Josephson junction (LJJ) devices ^{|. Another 
possibility that has been considered in the literature is 
that of nonlinear Klein-Gordon system where the on-site 
potential is ratchet-like [lq|. Notwithstanding, to our 
knowledge the case of spatial-symmetry breaking by in- 
homogeneities has not been studied in depth. One such 
study has been done by Carapella et al. who used an 
inhomogeneous magnetic field to create an effective in- 
homogeneous junction profile for fluxons to propagate 
[r?|. Recently an alternative to the generation of mo- 
tion for extended systems with a disorder in the chain 
has been proposed [18| . The novelty of the procedure is 
the design of a ratchet device from a lattice of point-like 
inhomogeneities. For this system, net motion arises from 
the interplay between disorder and nonlinearity of the 
nonlinear systems [l9|. 

In this paper we elaborate on the preliminary results 
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reported in |18| . Our aim is to carry out an in-depth 
analysis of the system, including a careful comparison to 
related point-like ratchets [2(J and an extension of our re- 
sults, originally obtained for the sine-Gordon (sG) model, 
to other nonlinear Klein-Gordon models such as the <f> 4 
equation. Additional motivation for this work arises from 
research on models of energy propagation along micro- 
tubule filaments inside the cells [21j . This application is 
specially interesting in view of the possible connection 
with the dynamics of transport in molecular motors in 
biological systems, with features similar to those of soli- 
tons as extended objects. In this context, the present 
work sheds light on the role played by the length scale 
competition between the point-like inhomogeneities (dis- 
order) and the size of kinks in the transport dynamics. 
For this purpose, we use the framework of collective co- 
ordinates (CC) in order to gain insight in the cause of the 
motion and the degrees of freedom that take part in it. 
Emphasis will be given to the kink width oscillations and 
their role in the transport properties: Indeed, in general, 
the width of the nonlinear topological excitations is cru- 
cial for the movement of these coherent excitations. The 
coupling between the translational and kink width de- 
grees is such that motion takes place 14]. In particular, 
in the present work we will see that the oscillation range 
of the kink width is determined by the interplay with 
the inhomogeneities. An additional, relevant issue is the 
analysis of the motion dynamics under thermal fluctua- 
tions. In ratchets, the noise is an important source of 
energy and, for some biological systems, it is regarded 
as the main cause of transport. Here we will consider 
the robustness of our rocking ratchet system under ther- 
mal fluctuations. In this case, activation of unidirectional 
motion was observed for a certain range of frequencies in 
the simulations as well as in the CC. 

In order to achieve the above mentioned goals, our pa- 
per is organized as follows: In Sec. II we formulate the 
basis for the ratchet device and explain the origin and 
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physical reasons of the rectification process. A discus- 
sion in the CC framework is devoted to the length scale 
competition between the inhomogeneities and the kink 
width, and its influence on the motion dynamics. In 
the same context, we establish an analogy with a sim- 
ple model used for describing unidirectional motion in 
molecular motors |20j, pointing out the relevant role of 
the kink width for the ratchet dynamics and its applica- 
tion in biological systems. Subsequently, an analysis of 
the efficiency in terms of the quantization of the trans- 
port is done for the relevant parameters of our ratchet 
system, including the interference effects among the in- 
homogeneities. An example is given for the <f> 4 model, 
not only for its known rich internal dynamics reinforced 
by the presence of an internal mode, but also for its po- 
tential application to macromolecules (e.g., in transfer of 
energy in microtubules |2l|). Next, in Sec. Ill, we analyze 
the kink dynamics subject to noise. In particular, a new 
phenomenon of activation of motion induced by noise is 
described. In order to explain this new behavior, dif- 
ferent collective coordinate approaches are implemented 
and thoroughly discussed. Finally, in the last Section we 
summarize the main contributions of our work and make 
a discussion about the perspectives on this line. We in- 
clude appendixes where we detail the CC approaches for 
one and two collective variables, extending the General- 
ized Traveling Wave Ansatz (GTWA) to the case where 
inhomogeneities and noise act together with damping and 
ac forces. 



II. RATCHET MODEL AND TRANSPORT 
A. Model 

Kink dynamics in the presence of inhomogeneities can 
exhibit different and interesting behaviors, depending on 
the interplay between those inhomogeneities and the non- 
linearity |22j,|23j among other factors. The generation of 
net motion using a lattice of point-like inhomogeneities is 
a good example of such non-trivial phenomenon |18| . Al- 
though in that previous work the problem was discussed 
for the sG model, it can be generalized in principle for 
any nonlinear Klein Gordon system. Therefore, aiming 
at that general viewpoint, for our analysis we formulate 
the model in a general way as follows: 

dU 

<j>tt + &<j>t - <Pxx + ut 1 + V(x) = Asm(ut + 5q), (1) 

0(j> 

where U(<f>) is the potential for the nonlinear Klein- 
Gordon equations and Asm(ujt + Sq) = f(t) is an ex- 
ternal ac force whose parameters A, uj and Sq represent 
the amplitude, frequency and phase of the periodic force, 
respectively. In particular, we will focus on the (/> 4 and sG 
models as specific examples; the corresponding nonlinear 
potentials are U{cf>) = \{<j) 2 -l) 2 and U(<p) = [l-cos(»] 
respectively. 



For our potential to exhibit ratchet-like phenomenon, 
we choose V(x) to be given by a periodically repeated 
unit cell, formed by an asymmetric array of delta func- 
tions (inhomogeneities). The unit cell configuration, of 
length L, is defined by three inhomogeneities, the first 
one located at the beginning of the cell, the second at 
distance a from the first, and the third at distance b from 
the second. The corresponding mathematical expression 
is 

V(x) = e [S(x — xi — nL) + S(x — X2 — nL) 

n 

+5(x -x 3 - nL)} (2) 

where the parameters satisfy the following constraints 
a, 6, c ~ Iq (static kink width in absence of inhomo- 
geneities); a, b < c with a ^ b, where L = a + b + c, 
a = X2 — Xi, b = X3 — X2 and c = L + x\ — X3, with 
xi < X2 < X3. For our study we have taken e > 0, where 
in the case of sG, specifically for LJJ t he p oint-like inho- 
mogeneities represent microshorts [24, 25] . However, the 
same scheme of arrays of point like inhomogeneities can 
be implemented for e < 0, where particularly for the sG 
model many works have been devoted 26] . The choice of 
three inhomogeneities is inspired by biological polymers 
like DNA where the existence of three bases per codon 
seems to be the ideal configuration for the occurrence of 
net transport [27|. In principle it is possible to obtain 
unidirectional motion by using an array whose configu- 
ration presents more than three inhomogeneities per unit 
cell if the distances between the delta functions are in the 
same length scale as the kink width (otherwise, different 
behaviors could arise like obtained e.g. in [23). How- 
ever, the inclusion of more inhomogeneities diminishes 
the efficiency of the transport as we will see below. 

B. Simulations of the model 

Contrary to the case of point particles, where motion 
through point-like inhomogeneities (delta functions) is 
physically meaningless, in our case we deal with kinks 
(extended objects) with a determined width. This is an 
intrinsic feature of these nonlinear excitations, and corre- 
spondingly the competition between their width and the 
distances among the inhomogeneities is crucial for the 
kink motion. The interference effects among the inhomo- 
geneities [28| create an effective potential for the motion 
of the kink center, the location of the inhomogeneities 
determining the direction of motion. For the particu- 
lar configuration of three inhomogeneities per unit cell, 
directional motion takes place only under the condition 
a =/= b. In the upper panel of Fig. ^ results of simulations 
of Eq. for the sG case with different values of a and b 
are depicted. Such a picture shows that our ratchet de- 
vice is an authentic rectifier. In addition, as in standard 
ratchet systems, the directional motion of the kink cen- 
ter takes place only for certain values of the amplitude of 
the ac force (see lower panel of Fig. a behavior that 
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FIG. 1: Simulations of Eq.Q for sG case: Position of kink 
center vs time: Upper panel: Different arrays with the same 
amplitude of the force A — 0.35: x\ = 0.5, xi = 1.8, 23 = 2.3 
(a > b) (solid line); xi = 0.5, £2 = 1.4, £3 = 2.3 (a = b) 
(dashed line); x\ — 0.5, X2 = 1, X3 = 2.3 (a < 6) (dash- 
dotted line). Lower panel: For different amplitudes of the ac 
force: A = 0.35 (solid line); A = 0.45 (dashed line); A = 0.50 
(dash-dotted line) with the array xi = 0.5, X2 = 1, X3 = 2.3. 
The other parameters used are uj = 0.05, e — 0.8, 8q = ir and 
period L — 4. 



is dependent on the ac force frequency. A more detailed 
picture of the dynamics of the mean velocity for the sG 
kink center as a function of the ac force amplitude for 
different frequencies can be found in |18| . 

We have restricted ourselves to the overdamped case 
by taking (3—1, where the inertial effects are small, re- 
ducing the generation and propagation of phonons, and 
for which the kink center moves on a tilted effective po- 
tential due to the external ac force. This regime pre- 
vents also the dependence on the initial conditions for 
the dynamics 29] . For the integration of Eq. JTJ we 
have used a Strauss- Vazquez numerical scheme |30] with 
free boundary conditions and spatial and temporal steps 
Aa; = 0.1 and At = 0.01, respectively. We have checked 
our results with different spatial steps Ax — 0.05 and 
Ax — 0.02. The simulations were done for the spatial in- 
terval [—30, 150] with inhomogeneities arranged periodi- 



cally according to our period in [0,120]. We have used 
the numerical representation for delta functions given by 



5(x - x c ) 



1/Ax, \x - x c \ < Ax/2 

0, otherwise. 



(3) 



Such representation is not unique but is one of the sim- 
plest forms for the numerical implementation p2l |2S| . 



C. Collective coordinate approach 

In order to understand the dynamics of the kink center 
different CC approaches have been implemented ^^|. In 
this previous work, it was shown that a 1-CC approach 
is not enough for the correct quantitative description of 
the results of the simulations even considering possible 
relativistic effects for high values of ac the force. How- 
ever, this approach allows us to understand the underly- 
ing physics of this ratchet system in the simplest way. On 
the other hand, satisfactory agreement with the simula- 
tions was found by an improved, 2-CC approach, which 
takes into account the kink width dynamics ^^] . At this 
point we will proceed directly with the discussion of the 
improved theory. The corresponding collective variable 
equations of two degrees of freedom for the Eq. (see 
Appendix B, Eqs. (|B2l) - llB22ll for D = 0) can be ex- 
pressed as 



X X 
M l j + PM l — - Mnl 



O'O' 



XI 



dU 
~dX 



if if), (4) 



aM l j +0aM o l o l - 
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-aM l - 



dU v 



dl 



— (5) 

dl ' (5) 



where X stands for the position of the kink center, I 
represents the kink width, the internal potential energy 
of the kink is 



jj-int 



1 



(6) 



and U(X, I) is an effective potential depending on the 
specific equation one is dealing with. In particular, for 
the sG case M = 8, Iq — 1, a = ir 2 /12, q = 2ir and the 
effective potential is given by 



U(X,l) = 2eJ2 



1 



+ - 



1 



cosli [(X-m -nL)/l] 
1 



cosh^ [{X -x 2 - nL) /I] cosh^ [{X - x 3 - nh) /l] 



■(7) 



As we can see from the previous equations the kink 
width dynamics is coupled to the motion of the center of 
the kink. Therefore, changes in the kink width directly 



FIG. 2: sG: Amplitude of the kink width versus kink center. 
Simulation (solid line); 2-CC Approach Eqs. 14151 (dashed 
line). The parameters are ui = 0.1, A = 0.44, for the array 
x\ = 0.5, X2 = 1, a?3 = 2.3, period L — 4 and e = 0.8. See 
text for a discussion of the loop. Inset: enlargement of the 
loop indicated by an arrow in the main figure. The motion of 
the kink center is indicated by numbered arrows. 



affect the translational motion. It is possible to observe, 
for instance, that decreasing the kink width decreases the 
effective ac force, making necessary to increase the am- 
plitude of the ac force in order to compensate for such 
effect. This is an important factor that explains in part 
the shift observed in the locations of the windows of mo- 
tion of the simulations with respect to those obtained 
from the 1-CC approach (see Fig. lb in 18]). Another 
relevant conclusion is the feedback between the effective 
potential landscape and the kink width, determined in 
turn by the potential. In this fashion, the 2-CC approach 
reflects the non-trivial interaction of the kink with the in- 
homogeneities, which is otherwise known to exhibit many 
counterintuitive phenomena [28}. 

To deepen our understanding of the dynamics, let us 
look into the oscillations of the kink width. As in the case 
of simulations we restrict ourselves to the overdamped 
case (taking (3 — 1). A picture of the kink width os- 
cillations versus kink center position is shown in Fig. [21 
As we can see, the agreement with the CC approach is 
excellent, confirming the validity of our predictions. It 
is particularly interesting the existence of loops, which 
arises as a consequence of the kink center motion back 
and forth into the wells of the effective potential (see Fig. 
|3J), i.e., describing half an oscillation before overcoming 
the barrier. 

Interestingly, another feature that stands out clearly 
is that the oscillations are around a value different from 
io = 1, the width of the unperturbed kink. Figure Ql 
shows that they take place around I ~ 0.8 and, further- 
more, that Iq is not even included in the range of os- 
cillations. This phenomenon is the result of the balance 
between two opposite forces. On the one hand, the inclu- 




4 6 8 10 12 

X 



FIG. 3: a) 3D plot of the effective potential Eq. Q. b) Frontal 
view rotated by a small angle. The array used is the same as 
in FigEl 



sion of inhomogeneities increases the potential energy of 
the system. This fact is reflected in the effective potential 
energy landscape Fig. [3J Such picture shows that when 
the kink width decreases, the potential energy decreases 
as well. Taking two points with the same value for X 
but with different kink widths I, for example M and O 
in Fig. Eb), we see this difference in potential energy, i.e., 
Em > Eo where Im > lo- Therefore, as the tendency of 
the system is to go to the minimum of the potential en- 
ergy, the kink width would decrease. On the other hand, 
the kink internal potential energy, Eq. ©, has a mini- 
mum in ioj and hence the energy increases when the kink 
width decreases (see Fig. (QJ for I < lo); notice that the 
first term of this equation accounts for a repulsive inter- 
action and the second is given by an attractive one. As 
a result of this balance, a new minimum will appear for 
the oscillations of the kink width. It is important to note 
that the difference AU mt of the internal potential energy 
for the kink width I = 0.7 with respect to the value Iq = 1 
(inset of Fig.QJ is of the same order as the energy differ- 
ence Em — Eo between the points mentioned before for 
the effective potential introduced by the inhomogeneities, 
in agreement with this discussion. 
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FIG. 4: Internal potential energy vs normalized kink width, 
Eq. © (sG case). The inset shows the part of curve where 
the oscillations of the kink width take place. 



D. Related point particle models. 

A problem closely related to our 2-CC approach, given 
by a point particle ratchet with two degrees of freedom, 
has been studied in [2(j. This model was designed for 
describing molecular motor dynamics consisting of two 
particles joined by a spring moving in a ratchet potential. 
The corresponding equations of motion are given by 



Ui 



U2 = —- 



dV(ux) dW(u 2 -ut) 
dV(u 2 ) dW(u 2 - ux) 



du 2 



du 2 



+ Asm(ut) + £ 1 (t), (8) 
+ Asin(ut) +&(*), (9) 



where Ui, u 2 represent the coordinates of the particles, 
V is a sawtooth potential and W is the internal potential 
energy. Here £j with £ = 1,2 are Gaussian white noises. 
Ignoring the noise terms and their influence on the net 
motion, we see that the change of variables X = + 
u 2 ) and I = u 2 ~ u\ casts the system in a similar shape 
as Eqs. (|4I5|I in the overdamped case, for which in good 
approximation the inertial terms could be neglected. In 
such a new context the variables X, I can be interpreted 
as the mass center and elongation (distance between the 
particles) respectively, and obviously they resemble to 
the variables mass center and width of the kink in our 
system. 

Notice that in both models we have an asymmetric po- 
tential. In our case it is given by Eq. J7|) , which is asym- 
metric at the CC level if the already mentioned condi- 
tions for the distances between the inhomogeneities are 
satisfied. In both systems, there are internal potential 
energies that characterize their elastic properties. In the 
model in |20| , the internal potential is expressed through 
a harmonic function (in the original variables): 



w(i) = h[i(t) 



(ii) 



i 



W{ui,u 2 ) = -k[{u 2 - ui) - l 



(10) 



where k is the elasticity constant. 

The links between the two models can be made more 
explicit by using a value for Iq close to the minimum 
around which the kink width oscillates in our simula- 
tions (cf. the discussion in the preceding subsection). 
It is important to point out that we have defined I as 
a kink width variable as given in the expression 1|B7|) 
of the appendix. However, what the quantity I actu- 
ally means is the distance at which the kink shape ap- 
proaches its asymptotic values, measured from the center. 
This means that I in our notation is half the "real kink 
width" . Consequently the ratio (real kink width) /(period 
of the effective potential) becomes 21q/L w 0.4 for which 
a very interesting dynamics for point particles dynamics 
has been reported in related 2-particle model |3lj |. 
This comparison between our model and that in |'20| al- 
lows to point out their main differences as well. It is par- 
ticularly important that in our framework, the internal 
energy can describe satisfactorily the repulsive interac- 
tion between real molecules where a van-der-Waals like- 
force prevents their overlap. This is very close to what 
occurs in molecular motors: if we take again the motion 
of kinesin as an example, this molecule has two dimcr 
heads that act as 'feet', allowing the molecule to 'walk' 
along a microtubule [32(. The repulsion would then ap- 
pear when the two dimer heads are too close. Such a 
repulsive interaction can not be naturally accounted for 
within the model of two particles. For solving this prob- 
lem the authors of [2(| resort to fix arbitrary values for Iq 
which in our case is not necessary. Note, however, that 
in spite of the technical differences between both models, 
phenomenologically they are very similar: both of them 
try to understand how the motion of molecular motors, 
which proceeds in steps accompanied by deformations (in 
the case of kinesin, when one step advances in front of 
the other) can arise. The common conclusion is that a 
point particle ratchet would not be a good model be- 
cause the second degree of freedom is needed to capture 
the whole mechanism of the motion. The advantage in 
our approach is that this second degree of freedom arises 
on its own, without a priori constructions, as an emer- 
gent property of the nonlinear excitation. Recent studies 
|3ll l33| show similar phenomena for the two degrees of 
freedom point particle ratchet of [20J when the ratchet is 
of flashing type. The close relationship of the model of 
[20| to ours suggests that nonlinear Klein-Gordon models 
can also exhibit rectification working as flashing ratchets, 
an issue we will address in future work. 



E. Length scales and quantization of transport 

It should be clear, from the results discussed so far, 
that in order to obtain a ratchet device for extended non- 
linear systems with topological nonlinear excitations, the 
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For the sG case we have the following expression 



FIG. 5: Normalized effective potential for the kink center 
within the CC approach, Eq. 1131 - 1141 . for two different ar- 
rays with e = 0.8. (a): xi = 0.5, X2 = 1, 003 = 2.3 and L — 4. 
(b): xi — 0.7, X2 = 1.4, x- A — 3.2, L — 5.6. In both panels sG 
(solid line); (j> 4 (dashed line). 



configuration of the inhomogeneities should be designed 
in such a way that the distance between the inhomo- 
geneities is of the order of the kink width. However, 
this picture is somewhat too simple, and as we will see 
below, another important factor to take into account is 
the existence of interference effects. Naively, one may 
try to design a similar ratchet system for the (f> 4 equa- 
tion. Considering only the kink width factor, it would 
seem that enlarging the sG array by a factor of V2 (the 
ratio between the kink widths in both models) similar 
phenomena would be observed. 

Let us make a more specific comparison between both 
models. To this end, we use the 1-CC framework in the 
nonrclativistic approach, where the equation of motion 
can be written as 



X + f3X = -^-^-sm(u J t + S ), (12) 



where u = U /Mq is the normalized effective potential. 



u(X) 



2e 



EE 



-{ cosh. 2 [(X 



iL)/l 



(13) 



with Iq = 1, and Mo = 8, whereas for the case of <fr 4 we 
have 



u(X) 



4M 



EE 



cosh [{X — Xi — nL)/lo] 



(14) 



with Iq = y/2, and Mq = 2\/2/3. The parameter q is the 
topological charge and it is given by q — 2n and q = 2 
for the systems sG and cf> 4 , respectively. The normalized 
effective potentials for two different arrays are depicted 
in Fig.|SJ Panel a) shows standard asymmetric potentials 
for ratchet systems obtained with an array that satisfies 
the conditions mentioned above for the location of the 
inhomogeneities in the sG case. However, in case b) the 
effective potentials obtained for an array approximately 
given by the multiplication of the factor y2 of the first 
one, shows a local minimum similar to an array of asym- 
metric double-well traps. This last one has been used as 
a device for generating motion of vortices in supercon- 
ductor materials j34[. 

According to the our previous arguments based on the 
kink width role, a similar picture is expected for the nor- 
malized effective potential of <fi 4 and sG if the arrays ver- 
ify the same length ratio as the full systems. Strikingly, 
Fig. [5] shows that the normalized effective potentials are 
almost the same but for the same array length. This ap- 
parent discrepancy can be explained if we take a detailed 
look at the potential given by Eqs. (|13I14|) for both cases 
(sG and cj) 4 ). It is clear from those expressions that, while 
in the case of 4> 4 we have a cosh 4 factor in the denomi- 
nator, sG has a cosh 2 factor. Therefore, the peaks and 
valleys in the effective potential for the cf> 4 system are 
much narrower than for sG, thus compensating for the 
increment in length. In addition, as in the sG model we 
will have dynamical changes of the effective potential due 
to the kink width variations, making more complicated 
the dynamics of motion. In any event, these effective 
potentials obtained in the simple approach highlight the 
importance of interference effects (see also |23|,|23]) and 
make it clear that the kink width is not the unique quan- 
tity to take into account. 

The consequences of choosing either the first or the 
second array for the kink dynamics are revealed in Fig. 
We have taken for the analysis the cj) 4 model with 
a relative low frequency for the ac force, for which the 
mean velocity as a function of the ac force amplitude 
shows a staircase structure. The range for the ampli- 
tude values was taken from the following rescaling ex- 
pression: q^A^/Alf = q sG A sG /M^ G . This rule is de- 
duced from the comparison between the 1-CC approaches 
for sG and <fi 4 models, considering the similitude of the 
normalized potentials discussed above. Fig. shows the 
dependence of the kink mean velocity as a function of 



7 




-0.005 
-0.01 

A 

§-0.015 
-a 
v 

-0.02 

-0.025 

-0.03 


-0.005 
-0.01 

A 

X -0-015 

-13 
V 

-0.02 
-0.025 
-0.03 



o.i 



0.2 

A 



0.3 



(d) 



0.1 



0.2 



0.3 



FIG. 6: Simulation results for 4> 4 : Mean velocity vs driving 
amplitude A for the frequency to = 0.015. (a) e = 0.8, xi = 0.5, 
x 2 = 1, x 3 = 2.3, L = 4(b) e = 0.6, xi = 0.5, a; 2 = 1, 
x 3 = 2.3, L = 4. (c) e = 0.8, an = 0.7, x 2 = 1.4, x s = 3.2, 
L = 5.6. (d) e = 0.6, an = 0.7, a; 2 = 1.4, x 3 = 3.2, L = 5.6. 
The lines connecting the points serve as guides for the eye. 



the amplitude for two different arrays and heights of the 
perturbations introduced by the inhomogeneities. As we 
can see, the motion is quantized as in standard ratchet 
systems 0, Q and is characterized by the existence of 
gaps for which the net motion is absent (i.e., pure os- 
cillating states). The mean velocity can be expressed 

as \(dX/dt)\ = \(V)\ = as usual [H, where the 

indices to, n S N quantize the motion. Using the expres- 
sion for | (V) | we can characterize the motion for each 
frequency and period of the array. Comparing the values 
obtained from the simulations with the results derived 
from the expression for | ( V) | with corresponding param- 
eters L and uj, we find that to and n can take the following 
values: For panel a), to = 1,2, 3, 4, 5 and n = 1; for panel 
b), to = 1,2,3,4 and n = 1; for panel c), to = 1,2 with 
n = 1 and n = 2, and for panel d), to = 1, 2 with n = 1 
and n = 2 . Although the mean velocity increases with 
the spatial period, the maximum value of the index to sig- 
nificantly decreases, leading to a global decreasing of the 
velocity. These results prove that the inclusion of more 
inhomogeneities per cycle, which obviously enhances the 
spatial period, is not a good option if we want to reach 
high velocities. Furthermore, a very low frequency would 
be required to obtain windows of motion. In the case of 
the dependence on the inhomogeneities height, the start- 
ing point of the stair-steps structure shows a shift to- 
wards greater amplitude of the ac force when increases 
the height, which is natural in order to overcome the bar- 
rier. Nevertheless, a higher speed is found, arising from 
a higher to and observable also as a broadening in the 
windows of motion. 



III. DYNAMICS UNDER NOISE 



This far, we have analyzed the ratchet-like behavior of 
our system in the deterministic case. However, it is clear 
that for our model to be more realistic, for instance, in 
the context of LJJ, the effect of the temperature has to be 
taken into account. The behavior of ratchet systems for 
nonzero temperature has been extensively studied both 
for point particles pi l3tj l37l l38j and for nonlinear 
extended systems [l(J u2, Ha]- As for the problem we 
are considering here, the fact that there has not been 
much effort on soliton ratchet-like phenomena induced 
by spatial inhomogeneities carries over to the stochas- 
tic effects. Therefore, it is important that we address 
this issue here. For the present work, we will focus on 
the robustness of our rocking ratchet under thermal fluc- 
tuations. Another relevant issue would be the possible 
activation, resonances or modifications of the transport 
features induced by noise, but this topic deserves a de- 
tailed analysis and will be the subject of future work. 
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FIG. 7: Mean kink velocity (dX/dt) vs driving amplitude A 
for different intensities of the noise, (a) ui — 0.05: circles 
D = 0., filled squares D = 0.005, squares D = 0.05; (b) 
uj — 0.1: circles D = 0., filled squares D — 0.005, squares 
D = 0.05. Lines serve as guides for the eye. The inset shows 
several realizations for the motion of the kink center with 
A = 0.43 , 5 = tv and D = 0.05. 



A. Full model 

For the sake of definiteness, we consider the sG model 
under the influence of a Gaussian white noise; the results 
for the </> 4 equation are similar. Introducing the effect of 
the temperature through the fluctuation-dissipation rela- 
tionship and considering the overdamped case as before, 
taking j3 = 1, we find the following equation: 

<i>tt +<k~ 4>xx + sin(0)[l + V{x)] = f(t) + ri(x, t), (15) 

with 

(Vfat)) = 0, 
(r)(x,t)T][x',t')) = DS{x- x')6(t-t'), (16) 

where f(t) = Asm(u)t+So) and the intensity of the noise 
D = 2k B T. 

For the numerical simulations of the full partial differ- 
ential equation as well as for the numerical solution of 



the collective variables (to be discussed in the next sub- 
section), we have used the Heun method with the Box- 
Muller- Wiener algorithm for generating Gaussian ran- 
dom numbers of mean zero and variance one |40j| . In 
Fig. 13 we show the behavior of the kink center dynam- 
ics under thermal fluctuations. Hereafter, we have set 
the array parameters to be X\ = 0.5, x-i = 1., 2:3 = 2.3, 
L = 4 and e = 0.8 for our study. The mean velocity was 
calculated using the expression in |8(, namely 

(V) -(*)-,*■ (*<" -*<"» , (!„ 

t— >oo t 

where the average is to be understood over many real- 
izations of the noise. From this figure we see that the 
steps of the deterministic case are now smoothed, a typi- 
cal feature for the dynamics under noise. It is important 
to realize that this smoothing affects the regions between 
the windows, which become minima of the mean velocity 
modulus |(V)| instead of gaps with null mean velocity 
(see Fig. UJi). This phenomenon is directly related to 
the strength of the noise, i.e, when the noise increases, 
the absolute value for the mean velocity decreases but 
simultaneously the connection between the windows be- 
comes more evident and the windows of motion become 
less pronounced. 

As in most other ratchet systems, in our model the 
stochastic fluctuations due to temperature assist the 
jumps of the kink center from one well to the next one, 
allowing in some cases jumps in the opposite direction of 
the rectification (see the inset graphic in Fig. |7Jd) which 
is not possible in the absence of noise. Accordingly, the 
thermal fluctuations affect the mechanism of rectification 
whereas, on the other hand, they yield the dynamics of 
the pure oscillating states of the kink center unstable (i.e., 
destabilizing the regions with locked directional motion 
for zero temperature). The joint action of both effects 
leads to the smoothing of the windows and the connec- 
tion of the deterministic gaps. For relative high tempera- 
tures the thermal kink energy is sufficient for overcoming 
the barrier of the effective potential, and the kink mo- 
tion is in practice diffusive, the influence of the barrier 
becoming negligible. For this reason the rectification of 
motion takes place only for not so large values of the 
noise intensity (see discussion in p4|). 

A remarkable feature we have observed in the simula- 
tions is shown in Fig. for frequency u> — 0.1, where 
new windows (absent in the deterministic case) appear. 
This scenario is very similar to the one reported in [l^ 
where a similar surprising and intriguing phenomenon 
was noted. There, the authors discussed that these new 
windows arose due to jumps of the fluxons between sta- 
ble and unstable pinned fixed points of the deterministic 
dynamics. Considering the interest for this stochastic 
phenomenon, we carried out a careful analysis of the cor- 
responding zone. To summarize this work, in Fig. [S] we 
plot the mean velocity as a function of the noise inten- 
sity for different values of the frequency, showing the ex- 
istence of an effective value for the intensity of the noise 
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FIG. 8: Mean kink velocity (dX/dt) vs intensity of noise D. 
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0.75. Inset shows one realization for the motion of the kink 
center for u> = 0.1, A = 0.70 , 5 = and D = 0.005. 



FIG. 9: CC approach: Mean kink velocity (dX/dt) vs driving 
amplitude A for different intensities of the noise and frequency 
uj = 0.1. Continuous line: D = 0, dotted line: D = 0.005, 
dashed line: D = 0.05. Inset: Mean kink velocity (dX/dt) vs 
intensity of noise D for A = 0.7625. 



for which a maximum absolute value of the mean velocity 
is obtained. The inset in Fig. [S] makes it clear that, as 
expected and suggested in |l2j, the mechanism of acti- 
vation occurs through jumps between multistable states 
(states of the kink center which in absence of noise are 
purely oscillating). Therefore, a higher velocity is ob- 
tained when the residence time in these multistates is 
reduced or, in other words, when the intervals between 
consecutive jumps decrease. Once again, this process of 
activation becomes more effective when the noise inten- 
sity increases, but above a certain value of the noise in- 
tensity the kink center starts to jump in the direction 
opposite to that of the rectification, leading to a global 
loss in efficiency. This explains the existence of an effec- 
tive value for the noise intensity for which the velocity 
reaches a maximum value. 

Another interesting characteristic observed in Fig.[S]is 
the dependence of the maximum mean velocity on the fre- 
quency. Specifically, for a frequency value slightly larger 
than u> = 0.1, the maximum velocity decreases, the peak 
moving towards greater values of the noise and the corre- 
sponding window of motion moving towards greater val- 
ues of the ac force. Accordingly, for a relatively large 
value of the frequency, above to = 0.11, the window of 
motion induced by noise disappears. On the other hand, 
for frequencies slightly smaller than u> = 0.1, a new win- 
dow in absence of noise is obtained. With all these re- 
sults, it is clearly established that the unidirectional mo- 
tion induced by noise occurs only for a narrow window 
of frequency values. 

We will show in the next subsection that this phe- 
nomenon seems to be a general feature, since at the CC 
level the system behaves very much like the dynamics of 
point particles. 



B. Collective Coordinates in presence of noise 

In order to understand the behavior observed in the 
previous section we resort again to the CC approach. As 
a first step, we take only into account one degree of free- 
dom. Although, as discussed above, this framework is 
inaccurate for describing quantitatively the kink motion 
on a lattice of inhomogeneities, it does help understand 
qualitatively most of the features observed in the sim- 
ulations, without unnecessary analytical complications. 
After some algebra (see Appendix A for details), with 
(3 = 1, we find the following expression for the kink cen- 
ter coordinate with noise: 

M X + M X = -^-qf(t) + v^DATo £(*) (18) 

with (£(t)) = 0, (£(*)£(«')) = S (t ~ For the sake 
of simplicity we have taken the nonrelativistic approach 
X 2 <C 1, for which the noise turns out to be simply ad- 
ditive. 

Figure |5| presents the results of the numerical integra- 
tion of Eq. (|18fl . Much as we did in the simulations, we 
calculate the mean velocity using Eq. I|17|) . taking up to 
500 realizations. From this plot two main features also 
observed in the simulations can be seen. First, smooth 
curves are obtained for the mean velocity as a function 
of the amplitude of the ac force, with values that de- 
crease when the noise is increased. Second, new windows 
appear, and inside them there is a value of the noise in- 
tensity for which the velocity reaches a maximum value 
(inset in Fig. 9). It is thus evident that, in spite of the 
quantitative differences with the simulations, this simple 
approach does predict correctly the qualitative behavior 
of the full system. 
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FIG. 10: CC Approach for two degrees of freedom (Eas. l|B210 - 
1B22I I of the appendix). Mean kink velocity (dX/dt) vs driv- 
ing amplitude A for different intensities of the noise, (a) 
(j — 0.05; continuous line: D = 0, dotted line: D = 0.005, 
dashed line: D = 0.05; (b) u) — 0.1; continuous line: D = 0, 
dotted line: D = 0.005, dashed line: D = 0.05. Inset 
shows Mean kink velocity (dX/dt) vs Intensity of Noise D 
for A = 0.72 with w = 0.1. 



In order to improve the results presented so far, we 
have extended the framework to two collective variables. 
By doing so (see Appendix B) we arrive at Eqs. (|B21|) - 
(|B22f) with two uncorrelated multiplicative white noises, 
meaning that the noises depend on the kink width dy- 
namics. The results for this improved approach are col- 
lected in Fig. ^3 Comparing with the simulations (Fig. 
0), we can observe the excellent agreement, with the loca- 
tions of the windows correctly predicted. As expected the 
curves are again smooth, a feature correctly accounted 
for already in the 1-CC framework. For the frequency 
u> = 0.1, a new window is predicted, whose location is also 
in very good agreement with its corresponding window 
in the simulations. These results confirm the importance 
of considering the kink width dynamics in the framework 
of the collective coordinates in order to achieve correct 



quantitative results as compared to the simulations. 



IV. CONCLUSIONS 

We have studied the dynamics of solitons in a lattice 
of point-like inhomogeneities in two different nonlinear 
systems, the sG and the (/> 4 models, as good represen- 
tatives of an entire class of nonlinear models. We have 
designed a periodic lattice made of an asymmetric unit 
cell, where, depending on the lattice parameters, ratchet- 
like motion of solitons (rectification) is observed. This 
is one of the few examples proposed so far for soliton 
ratchets based on spatial inhomogeneity. Building on a 
preliminary work reported in [Tsj . we have analyzed in 
full detail all the characteristics and features of interest 
of our model. In order to understand the observed phe- 
nomenology, we have developed two CC approaches: one 
that takes into account only the motion of the kink cen- 
ter, and another one that includes the kink width as a sec- 
ond degree of freedom. We have seen that whereas both 
techniques give good qualitative results and allow for a 
general understanding the phenomenon, only when the 
kink width enters in the CC description the results are 
quantitatively correct. Interestingly, the CC approach 
allows to show that rectification takes place only when 
the unit cell gives rise to an effective potential of the 
same length scale as the kink width. This fact, along 
with the relevant role played by the width dynamical be- 
havior, makes us conclude that one essential ingredient 
for observing ratchet-like behavior in this system is the 
existence of length scale competition ^| . 

As a necessary complement of the deterministic study, 
we have considered the influence of thermal noise on the 
behavior of our model systems. As in the determinis- 
tic case, the CC technique yields very good results. Of 
particular relevance is the appearance, at the 2-CC level, 
of multiplicative noise in the CC equations, coupled to 
the kink width. This is yet another hint to the crucial 
role played by this degree of freedom in the dynamics 
of the system. Another important feature is the motion 
activation induced by noise, that was observed in the sim- 
ulations and later corroborated in the CC framework. 

In a wider context, we believe that the CC approach 
presented here, and the conclusions about the role of the 
kink width, are of a much more general applicability. 
This is the case, for instance, with the anomalous res- 
onances observed in nonlinear Klein-Gordon models [42| . 
which can only be explained by considering the width 
degree of freedom. Specifically within the context of rec- 
tification, we have previously shown that the existence of 
ratchet behavior induced by pure asymmetric temporal 
driving in a homogeneous system, recently demonstrated 
in experiments in L J J , is due to the coupling between 
the translation and the width oscillations [l4| . Based on 
all these experiences, we believe that the implementation 
of the 2-CC approach can be very useful in the case of 
spatially correlated noise (IH Eij . where the 1-CC ap- 
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proach can not explain the length scale competition be- 
tween the correlation length and the size of the kink, 
first suggested in 0|. The new formulation would be 
also necessary in using the CC framework for describing 
the transport of proteins assisted by a thermal bath pro- 
voked by ATP molecules hydrolization, where the width 
is an important quantity to take into account. On the 
other hand, it is becoming more and more evident the 
crucial contribution of internal degrees of freedom in the 
functioning of molecular motors [32j. In this respect, our 
results suggests that models including in a natural way 
this degree of freedom, like in our case, can be the proper 
descriptions of those phenomena. Size dependent rectifi- 
cation has been also reported in colloids [43 and rectifi- 
cation through clustering has been observed in granular 
gases |46|. which further reinforce our conclusion that 
the deformation or internal degrees of freedom must be 
an ingredient of good theoretical description of ratchet 
phenomena beyond the point particle scenario. 

Finally, we want to stress that our design of a soliton 
system with ratchet behavior is a very simple one, that 
can be implemented in actual experiments and devices 
such as LJJ, for instance. Another field where the cur- 
rent state-of-the-art allows to apply this result is that of 
engineered molecular motors, demonstrated, e.g., in |47| . 
In this case, our proposal may apply to the design of 
biomolecular devices with medical applications. Gener- 
ally speaking, this type of approach to rectification can be 
of interest for applications in which it is needed to have a 
tunable rectifier tailored for a specific regime. Of course, 
experimental verification of our predictions is needed to 
ascertain the accuracy of our results. We hope that this 
work stimulates experiments in this direction. Indeed 
the research reported here, opens new perspectives in the 
design of ratchet devices for more complicated extended 
nonlinear systems, such as general coupled chains Eg- 
Of particular interest in this class are stacked LJJ |49| . 
although there are many other systems with potential 
applications in different areas. Work along this line is in 
progress. 
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where ip = ip, f it) = Asm(ujt + So), D = 2{ik B T and 
H is the Hamiltonian corresponding to the unperturbed 
form of Eq. given by 



h = J dx{-^ + -4>l + u{4>)}. 



(A4) 



As starting point we assume that the solution has the 
form 

4>(x,t)=<t> K [x-X(t),X], (A5) 
and therefore by definition of ip we have that 



rl>(x,t)=il> K [x-X(t),X{t),jt\. 
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The index K refers to the kink shape, but in the fol- 
lowing we will omit it for simplicity. 

Following the procedure in |50| . inserting p 7 up into 
Eas. (|Alf >- (|A2|) we get the expressions 
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APPENDIX A: COLLECTIVE COORDINATES. 
FIRST APPROACH 



dip ■ dip ■■ dip ••• 

-^x + —^x + —^x 

OX dX dX 



dU 



Sep f \dx ax 



V(x)+f(t)+r)(x,t). (A8) 



Multiplying Eq.^XTJ by and Eq.ljUJ by J^, and 

then subtracting both expressions and integrating we ar- 
rive at the following equation 



In order to apply the GTWA first proposed in |5C 
rewrite Eq.Q adding Gaussian white noise 



Ax + MX = -/3CiX-/3C 2 X+F ac +F stat +F mn +F st , 

(A9) 
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whose values for the coefficients and forces are given 



Hence the correlation function for F st can be written 



by 



r, = / ,i, [ ^ ) F inh 



— oo 
oo 



Co = / dx 
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dx^Vix)^ 

dt> ( >dx' 



dcf> 
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/OO 
dx T](x,t) 
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S<p dX Sip OX J 



where £ represents the energy of the system, TL is the 
Hamiltonian density of Eg . (|A4|) and F stat is the static 
force due to the external field, equal to zero for the above 
Hamiltonian. 

Next we consider the sG potential for the system 
Eas. HAl|) - (jA2|) for which we assume as solution the ansatz 

cp(x,t) = cj) {0) {"{(x-X(t))} =4arctan(exp[7(a;-X(t))]), 

(A10) 

where 4>^ = 4 arctan (exp \(x — Xq)/Iq]) is the static 
kink solution of the sG system, centered in Xq and of 



width Iq. Here 7 = 1/v 1 — X 2 where we have put lo = 1 
for the sG case. 

Considering the previous statement for the static force 
and taking into account V(x) from Eq.J2J), we obtain 



A = 0, F ac = -qf(t), 

), 

dU 



M = j 3 M , F stat = 



Ci = 7 M , F 



ink 



dX' 



c 2 



0, 



where Mq = 8 is the kink mass, q = 2tt is the topolog- 
ical charge and U(X, X) given by 



U(X,X) = 2eJ2J2 



1 



7? t~l C0Sh \?l( X ~ X i~ nL )\ 



(All) 



is the effective potential. In the non-relativistic limit 
X 2 < 1, U(X,X) ~ [/(X). 

A representation for the stochastic force F st can be 
obtained from the calculation of the variance. In the 
case of additive noise it is allowed to make the following 
assumption 



,d<j>W{x,t) d<f>W(x',t') 
v dX dX 

d(f>(°\x,t) d<p(°\x',t') 



dX 



OX 



(ri(x, t) V (x',t')).(A12) 
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(A13) 



for which, taking into account the expression l|A2|> . after 
some algebra we get 



(F 8t (t)F st {t')) = 2(3k B TjM 5(t - t% 



(A14) 



i.e., F st (t) is a white noise with kink diffusion constant 
D K = 7 M D. 

As a consequence we obtain a non-additive noise due 
to the factor j(X), i.e, we have a problem with multi- 
plicative noise. 

Then the equation of motion i|A9() can be rewritten as 



YM X + (3jM X 



-»/«)- g 



^ £(*) (A15) 



with (£(*)) = 0, (f(t)£(f)) = *(*-*')• The Eq.JSH 
in absence of inhomogeneities and noise agrees with the 
results presented in 51] . The other r.h.s. terms that ap- 
pear in (| A 1511 are in correspondence with those already 
obtained in |23t 152^ in the presence of impurities (non- 
relativistic approach) and Gaussian white noise, respec- 
tively. The procedure used here is equivalent to the so- 
called adiabatic approach by using modified conservation 
laws El. 



APPENDIX B: COLLECTIVE COORDINATES. 
SECOND APPROACH 



In order to get the CC equations we follow a similar 
procedure as in the previous section but this time we 
propose a solution with the form 



(x,t) = <p[x - X(t),l(t)], (Bl) 



iP(x,t) = ip[x - X(t),l(t),X(t),l] (B2) 

with ip = (p, which considers the kink width as a new 
collective variable (see e.g. [42]). 

Inserting Eas. (|Bl|) and l|B2|) in our system Eas. i|Al|) - 

(|A2|) and then multiplying the first equation by and 

the second one by Jh^- ; subtracting both expression and 
integrating we arrive at the following equation 
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On the other hand, Ea. (|B6|) is transformed into 
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where 7i is the Hamiltonian density of Eg. l|A4|) for 
which, as was seen before, a null value for F stat is ob- 
tained. 

Repeating the same procedure, but now with and 
we get the expression 



dx[ip, <p]X - 
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Following Rice [53 for the particular case of sG 



(p(x,t) = (p {0) [x~X{t),Ht)] = 4arctan | exp 
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(B5) where a = 7r 2 /12, Al = 8, Iq = 1 and 
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As in the previous section we use the variances of the 
stochastic forces in order to obtain approximate expres- 
sions for it. Taking the assumption given by the expres- 
sion 1JA12I) we find for (|B11|) the correlation function 
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In what follows similar expressions to the Ea. l(AT2| 
valid for additive noise are used in order to calculate other 
correlation functions like 
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, ,d<p^(x,t) d<p^(x',t') . . . . . ... 



— oo ./ — OO 



9Z 



DS(t-t') j dx(?^pj =DS(t-t')jaM , 



(B19) 



and 

(F st {t)K st (t')) 



dxdx' v y > v ^ 1 (vixMx 1 ,t')) 

oo J — oo C.A 



dl 



D5{t - t') 



DO- 



dx dx di =°- 



(B20) 
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From the latter correlation for the stochastic forces we 
see that these are not cross-correlated. 

Finally, collecting all the previous results we can 
rewrite Eas. (|B8f) . (|B13|) as follows 



MqIoj+/3Mq1 j 



, XI 



DMok 



I 



&(*)> (B21) 



. . . X 2 _ dU 



K int (l,l,X) 



DuMqIq 



I 



f 2 (i) (B22) 



with (&(t)) = o, = M(* - 0. fOT i = 

1,2. 

A feature of particular interest in these new equations is 
the presence of stochastic forces which are multiplicative 
white noises dependent on the kink width variable. 



The method described here using the technique of pro- 
jection is equivalent to the variational calculations of the 
momentum and the energy of the system for perturbed 
nonlinear Klein-Gordon systems of the form of Eqs.jAdJ 
and ()A2|) and with a Hamiltonian of the form of Eq. (|A4l) 
(see |42j for details). Another procedure and derivation 
has been recently presented in |54| . 
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